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Abstract. In this paper we prove uniform a priori estimates for transmission problems 
with constant coefficients on two subdomains, with a special emphasis for the case when 
the ratio between these coefficients is large. In the most part of the work, the interface 
between the two subdomains is supposed to be Lipschitz. We first study a scalar transmis- 
sion problem which is handled through a converging asymptotic series. Then we derive 
uniform a priori estimates for Maxwell transmission problem set on a domain made up of 
a dielectric and a highly conducting material. The technique is based on an appropriate 
decomposition of the electric field, whose gradient part is estimated thanks to the first part. 
As an application, we develop an argument for the convergence of an asymptotic expansion 
as the conductivity tends to infinity. 



1. Introduction 

The goal of our work is to derive uniform a priori estimates for transmission problems 
in media presenting high contrast in their material properties. We investigate in particular 
the heat transfer equation 

(1.1) div a grad (yj = / 

and the Maxwell equations given by Faraday's and Ampere's laws 

(1.2) curl E — zwyUoH = and curlH + {icueo — a)E = j . 

Here, a represents the heat conductivity and a the electrical conductivity. We assume 
that these equations are set in a domain Q made up of two subdomains fi+ and f2_ in 
which the coefficients a and a take two different values (a+, a_) and (cr+, cr„), respectively. 
These equations are complemented by suitable boundary conditions. Our interest is their 
solvability together with uniform energy or regularity estimates, namely 

• when the ratio |a_|/|a+| tends to infinity in the case of eq. (1.1) 

• when a+ = (insulating or dielectric material) and cr_ = a tends to infinity (highly 
conducting material) in the case of eq. (1.2). 

We address different, though connected, issues for these two problems, namely the issue 
of uniform piecewise regularity in Sobolev norms for solutions Lp of equation (1.1), and the 
issue of uniform estimates for the electromagnetic field (E, H) solution of system (1.2). 
None of these questions have obvious answers, all the more since we do not assume that 
the interface S between fi+ and Vt^ is smooth. 

Our whole analysis is valid under the only following assumption on the interface 

(1.3) S is a bounded Lipschitz surface in M^. 
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In the Maxwell case, similar estimates as ours are obtained in [7], but under a stronger 
regularity assumption on S. Our approach differs, being based on a decomposition of the 
electric field given in [2]. The gradient part of the decomposition is handled through the 
uniform regularity estimates proved for equation (1.1). 

The paper is organized as follows. In section 2 we introduce the notations and give the 
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main results. In section 3 we prove uniform piecewise Ha estimates for solutions of the 
scalar interface problem (1.1) with exterior Dirichlet or Neumann boundary conditions. 
In section 4 we prove uniform estimates for the electromagnetic field (E, H) solution of 
the Maxwell system (1.2) when the conductivity a of the conducting part is high. We 
conclude our paper in section 5 by an application of the previous uniform estimates to the 
convergence study of an asymptotic expansion as the conductivity tends to infinity. 



Let f2 be a smooth bounded simply connected domain in R'^ with boundary dVt, and 
r2_ CC be a Lipschitz connected subdomain of f2, with boundary S. We denote by 
the complementary of Vt^ in Vt, see Figure 1 . 



We denote by tp'^ (ip ) the restriction of any function (p to 17+ (f^-)- 

2. 1. Scalar problem. We consider both Dirichlet and Neumann external boundary condi- 
tions associated with equation (1.1) and introduce the functional spaces suitable for their 
variational formulation: Vd = HQ(f2) for Dirichlet and Vm = {y? G H^(f2) | J^(pdx = 0} 
for Neumann. For any given function a = (a+, a_) determined by the two constants a± on 
Q,± and in either case (V = Vn or V = Vd) the variational problem is: Find if E V such 
that 



2. Notations and main results 




Figure 1 . The domain f2 and its subdomains 17+ and f2 




where the right-hand side (/, g) satisfies the regularity assumption 
(2.2) / G L\n) and g G L\J:) 
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and the extra compatibility conditions 

(2.3) / / dx = and / ^ ds = if y = and gds = if V = Vd. 

Our main result in the scalar case is the following piecewise H^/^ a priori estimate, uniform 
with respect to the ratio p := a_(a_|_)^^. It applies both to Dirichlet and Neumann boundary 
conditions. 

Theorem 2.1. Let us assume that a+ ^ 0. There exist a constant po > independent ofa+ 
such that for all a- & {z ^ C \ \z\ ^ po|a^+|}. the problem (2.1) with data (/, g) satisfying 
(2.2)-(2.3) has a unique solution Lp G V, which moreover is piecewise H^/^ and satisfies 
the uniform estimate 

(2.4) y'-lln^ + y-\\l.n^ ^ C,,{\a+\-'\\f\\o,n + I^Ho.e) 
with a constant Cp^, > 0, independent of a^, a_, f, and g. 

This statement is proved in the next section using an asymptotic expansion for ip with 
respect to the powers of p^^ = a+(a_)^^. The estimate (2.4) will be a consequence of the 
convergence of this series in the piecewise H'^/^-norm. The dependence of po and Cp^^ on 
the overall configuration is discussed in Remark 3.1 after the proof. 

Remark 2.1. The estimate (2.4) is uniform for fixed a+ when |a„| tends to infinity. The 
roles of a+ and a_ can be exchanged and an estimate similar to (2.4) proved. In fact, there 
holds a more precise estimate where a+ and a_ play symmetric roles, see Proposition 3.3. 

Remark 2.2. In the Neumann case, the compatibility conditions (2.3) are necessary for the 
right hand side of problem (2.1) to be compatible for all values of (a_, a+), because of the 
factor (a_|. — a„) in front of the integral on S. If this factor is replaced by 1, then, under 
the weaker conditions 

(2.5) - fdx+ gds = if V = Vn and nothing if 1^ = Vd , 

Jn Jt, 

the problem (2.1) is still solvable for p large enough, see Proposition 3.4. 

Remark 2.3. The assumption that S is Lipschitz is necessary: There exist non-Lipschitz 
interfaces such that estimate (2.4) does not hold. In two dimensions of space, such an 
example is provided by the checkerboard configuration (Figure 2), cf. [5, Theorem 8.1]. 




Figure 2. A non-Lipschitz interface S 
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Remark 2.4. If the Lipschitz interface S is polyhedral, there hold uniform piecewise H** 
estimates for any exponent s, | < s ^ s^, with some | < ss ^ 2, cf. [14, Ch. 1, §1.5]: 

(2-6) llv'^l|s,n+ + y'\\s,n- ^ Cp^^s{\a+\~^\\f\\s-2,n + II^IL-|,s)- 

This is a consequence of our proof (see Remark 3.1) in combination with elliptic estimates 
in polyhedral domain, cf. [6]. In particular, if S contains some edge, then < 2. 

Remark 2.5. If the interface S is smooth, there hold uniform piecewise H'' estimates for 
any s ^ 2: 

(2.7) ii^^ikf^+ + ib"ikj7_ ^ (^po,4k+ri/L-2,n + ii^iL-|,s)- 

This is again a consequence of our proof in combination with standard elliptic estimates, 
cf. [1] for instance. 

2.2. Maxwell problem. We consider two types of boundary conditions to complement 
the Maxwell harmonic equations (1.2) on Either the perfectly insulating conditions 

(2.8a) E ■ n = and H x n = on dn , 

where n denotes the outer normal vector, or the perfectly conducting conditions 

(2.8b) E X n = and H ■ n = on dQ . 

In both cases, for the conductivity a = (0, a), we can prove uniform a priori estimates for 
the electromagnetic field as cr — t- oo provided the following condition on limit problems in 
the dielectric part ^2+ is valid: 

Hypothesis 2.2. The angular frequency u is not an eigenfrequency of the problem 

{curl E — zwyUoH = and curl H + zweoE = in VLj^ 
Exn = OanJH-n = onS 
(2.8a) or (2.8b) on dVt. 

Our main result for Maxwell equations is the following a priori estimate, uniform as 
(7 — 7- oo. The right hand side j is chosen to belong in H(div, Vt) or Ho(div, Vt) where: 

(2.10a) H(div,fi) = {u G L2(f^) I <livueV{VL)} with V{VL) = l}{VLf 
(2.10b) Ho(div, []) = {u G H(div, Vl) \ u ■ n = on dVL}. 

Theorem 2.3. Under Hypothesis 2.2, there are constants (Jq and C > 0, such that for 
all a ^ (To, the Maxwell problem (1.2) with boundary condition (2.8a) and data j G 
Ho(div, VL) has a unique solution (E, H) in L^(f2)^, which satisfies: 

(2-11) ||E||o,f7 + ||H||o,n + ||E||o,n_ ^ C'||j||H(div,f^)• 

A similar result holds for boundary conditions (2.8b) and data j G H(div, Vt). 

This theorem is proved in section 4. It is based in particular on a decomposition of the 
electric field into a regular field w in H^(r2) and a gradient field Vv? for which Theorem 
2. 1 will be used. 
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3. Proof of uniform scalar regularity estimates 



3.1. The problem. Under the assumptions of Theorem 2.1, we normalize the equations 
by dividing by a+. Denoting the quotient a_(a+)~^ by p, and still denoting by / the new 
right hand side (i.e., the old one divided by a+), we can write problem (2.1) in the form of 
the following transmission problem 

pAv9~ = /~ in f2_ 

= /+ in fi+ 

(3.1) <J v'" = ip+ on S 

9n¥?+ - pd„ip- = (l- p)g on S 

(b.c.) on 

where d„ denotes the normal derivative (inner for outer for (]+). The external boundary 
conditions (b.c.) are either Neumann or Dirichlet conditions. 

Our method of proof for Theorem 2. 1 consists in the determination of a series expansion 
in powers of for (p solution of (3.1): We are looking for solutions in the form of power 
series 



(3.2) 



n=0 

oo 

n=0 



in f2_ 



in Q. 



Since the expansions are different according to the boundary conditions, we treat first the 
Neumann case in subsection 3.2 and the Dirichlet case in subsection 3.3. We prove com- 
plementary results in subsection 3.4. 

3.2. Neumann external b.c. Inserting the ansatz (3.2) in the system (3.1), we get the 
following families of problems, coupled by their conditions on S: 



(3.3) 




= 


in 




= 9 


on 




and 












{ A^o 


= r 


in 




(3.4) 




= 


on 


s 






= 


on 


dQ 


and for k e W (here Sj, 


is the Kronecker symbol) 






(3.5) 


f A(p^ = 
I = 


Sir 

-Sl9 + d„ 




in 
on 


and 
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(3.6) 




= "Pk 


on 
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on 


dQ 
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Thus we alternate the solution of a Neumann problem in and a mixed Dirichlet- 
Neumann problem in $7+. Since we have assumed that S is a Lipschitz surface, we have a 
precise and optimal functional framework to describe these operators and their inverse. 
We need the notation (s < 2) 

(3.7) W{U, A) := {if G B'{U) \ Aip G V{U)}. 

Using [9, cor. 5.7] (see also [3] for a similar context) with the fact that S is Lipschitz, we 
obtain the following equalities between spaces on : 

(3.8a) {if d„ip\^eL\i:)} = Hi(l]_,A) 

(3.8b) {(^ G H^(fi_,A), (^1^ G H^(S)} = Hi(l]_,A) 

and on f2+ 

(3.9a) {^e}l\n+,A), ip\^eli\J:), d„ip\g^eL\dn)} = Hi(fi+,A) 
(3.9b) {^eE\Q+,A), d„^\^eL\j:), if\^^eE\dQ)} = Hi(fi+,A). 

As a consequence of the previous equalities, the following definitions for the resolvent 
operators and fH+ make sense and define bounded operators: 

• is the resolvent of the Neumann problem on 

(3.10) : |(F,G) G L2(fi_) X L2(S) I j F dx + j G ds = 0^ 

— > |$GHt(fi_,A)| J <l>dx = o} 

where $ = G) satisfies A$ = F in fi„ and d^^ = G on E, and 

• 9^+ is the resolvent of the Dirichlet-Neumann problem on 17+ 

(3.11) : |(F,i7) G L2(fi+) X Hi(S)} Hi(fi+,A) 

where $ = H) satisfies A$ = F in $ = i7 on S and = on dn. 

A further consequence of equalities (3.8)-(3.9) is that the following trace operators make 
sense and define bounded operators: 

• 1° is the Dirichlet trace on S from inside (]_:(/?" i-;- |s and is bounded 

(3.12) T°:Hi(0_,A) ^ H^(S) 

• is the Neumann trace on S from inside 1]+: h-> dn^'^lj: and is bounded 

(3.13) 1^:Hi(fi+,A) ^ L\E) 

Note that none of the two operators T° or %\_ is bounded if Hi A) is replaced by the 
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larger space Ha 
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3.2.1. Discussion of elementary problems. The Neumann problem (3.3) admits the solu- 
tion ip^ = y{_{0,g), since the compatibility condition J^g ds = holds by assumption. 
We have the estimate 



(3-14) Ibn II 3, , ^ C„||0 



where C_ is the operator norm of the operator . 

Thus belongs to H2 A). Hence its trace f/^o = V^o Is belongs to H^(S). Next, 
problem (3.4) admits the solution (p^ = and if^ belongs to Ht(f2+, A) 

with the estimates 

(3.15) llv'o llHi(n„A) ^ + 11^0 lli,s)- 

Here C+ is the operator norm of the operator 

Then we continue with problems (3.5) and (3.6) in a similar way, the only point to 
discuss being the compatibility condition in the Neumann problem (3.5). 

Lemma 3.1. The Neumann problem (3.5) is compatible. 

Proof. For k = 1, we must show that 

(3.16) / f-dx + [i-g + d„y^+) ds = 0. 

According to (3.4), Ai/Pq = /+ in Vt^ and dnip^ = on dVt. Integrating by parts, we get 



/ dn^-^ ds= /+dx . 



Thus, we deduce from hypothesis (2.3) the compatibility condition (3.16). 
For k ^ 2, let us assume that the term </?^_i was built. We must show that 

(3.17) / d„ipt_,ds = 0. 

Jt. 

According to (3.6), Av9^_^ = in and d„ip~l_^ = on dVl. Integrating by parts in 
we get (3.17). □ 

Consequently, for all A; ^ 1, the Neumann problem (3.5) admits the solution := 
—5\g + T^v9^_J. Then (/)^ belongs to H2 A) with the estimates 

(3.18) Ibfc llHi(n.,A) ^ C'_(5m|r||o,C_ + |k||o,E) + ||5nV^til|o,s). 

Finally problem (3.6) admits a unique solution Lp^ G Hi(r2+, A) with the estimate 



(3-19) ll^^llHi(.,,A)^^+ll^^"lli. 



s 



In (3.18) and (3.19), the constants C_ and C_|_ are the same as in (3.14) and (3.15). 
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3.2.2. Uniform estimates. Let C° and > be the operator norms of the Dirichlet and 
Neumann traces T° and l^, respectively cf. (3.12) and (3.13). We set a = C+C^C-C\ 
with the constants C+ in (3.15) and C_ in (3.14). According to (3.18) and (3.19), we see 
by an induction on ri G N* that 



(3.20) 



II 3 ^ a" ^||(j9i II 3 

l^^"llHl(n„A)^^-^°""-1v^rllHl(._,^)- 

-1, 



Let po > such that Pq a < 1. Then, for all p G C such that |p| ^ po, the series of general 
terms p~"</?~ et p~"</?J converge respectively in Hi(f2_, A) and Hi(f2+, A). We denote 
by (/)(p) the sum of these series. Moreover, normal convergence is geometric with common 
ratio |p~^|a, bounded by Po^a. Hence 

(3.21) 



According to (3.18) for k = 1, and (3.14)-(3.15) for k = 0, 

(3.22) |brllHl(n_,A)^^-(ll/"llo^f^- + ll^llo.^ + ^'^+(ll/''llo'^+ + ^°^-ll^llo^^))- 
With (3.14), (3.15), (3.21) and (3.22), we deduce the uniform estimate for |p| ^ po 

(3.23) ll<p)ll|,n^ + llv^(p)ll|,n_ ^ C{po) {Wfhn + ||^^||o,s). 

3.2.3. Proof of Theorem 2.1 in the Neumann case. By construction, </?(p) is solution of the 
problem (3.1). Hence, (p(^p) G H^(r2). Setting 



"P'ip) ■= "Pip) 



we obtain a solution cp = (/?'^^^ of the variational problem (2.1), which moreover satisfies 
estimates (3.23), hence estimates (2.4). 

It remains to prove that the solution of the variational problem (2.1) is unique when 
IpI ^ Po. Let (p^ eV he solution of problem (2.1) for such a p and for f = 0, g = 0: 



(3.24) WtpeV, / a+V(p+ • dx + / a_V<^; ■ dx = 0. 

Jn+ Jq- 

On the other hand, the power series construction yields a solution ^/^^, of problem (2.1) with 
a± instead of a± and with / = = (note that these data satisfy assumption (2.3)): 



(3.25) G 1/, / a+V^Z ■ dx + / a_VC " V<^" dx = - / ^dx, 

Jn+ Jn^ Ji 

Taking the conjugate of (3.25) for ip = ^p^, and (3.24) for ip = ijj^,, we find 

ip^ Tpldy. = 0. 

Hence the uniqueness, which concludes the proof of Theorem 2. 1. 
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Remark 3.1. From the above proof, we can see that the constants po and Cp,, depend only 
on the four operator norms C°, (trace operators T° (3.12) and %\ (3.13)), C_ and C+ 
(resolvent operators (3.10) and (3.11)). The extension of the estimates (2.4) to 
different sets of Sobolev indices, cf. Remarks 2.4 and 2.5, depends on the boundedness of 
the four operators (the orthogonality conditions are understood for the last two ones): 

1° : WiVt^) — ^ ff"5(S) , %\ : (fi+) ff^i(S) , 

: ff-2(fi„) X ff^i(S) ^ (fi„) , fH+ : ff^^^^^) ^ ff-^(S) ^ . 
In particular, none of them is bounded for s = |, so we cannot set s = | in estimate (2.6). 

3.3. Dirichlet external b.c. Here V = HK^l). When we consider the boundary condition 
= on d^l in problem (3.1), a similar construction can be done. However, we need 
a special care to treat the compatibility conditions in Vl^. Starting from the same Ansatz 
(3.2), we get 



(3.26) 
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in 
on 


S 


and 


















u 


in 




(3.27) 1 








on 


s 




< ¥>o 







on 


dQ 


and for /c = 1, 2, ... 












(3.28) 1 f^'^i 


-SI 


9 + 




in 
on 


and 














■ A^+ 







in 




(3.29) 1 








on 


s 











on 


dQ. 



3.3.1. Discussion of elementary problems. Let <^q G H2(f2_, A) be the solution of the 
Neumann problem (3.26) under the condition <^q dx = 0. Here we still keep a constant 

to be adjusted; we call it cq. Then (/Jq = <^q + cq will be determined once cq is fixed and 
(3.27) will give a unique v?o ^ Ha A). 

We consider now (3.28) for A; = 1, which is a Neumann problem with the compatibility 
condition 

/ /-dx+ / {-g + d„^^)ds = Q 

and since J^gds = 0, it reads 

(3.30) / d„ip+ds = ~ f f-dx. 
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But now we choose lPq=(pq+ co^Jj where 

= /+ in Q+ 
(3.31) { = <^o on S 




on on. 



and 




= in n+ 

(3.32) { =1 on S 

= on on. 

with 

(3.33) """^'iL "^'^ L ^ / L ^"^ 

Clearly (^q G H5(r24,, A) is uniquely determined by (3.31), by (3.32), and cq by (3.33) 
since J^dntp^is ^ (principle of maximum). 

Thus we have completely determined (/^o = (pQ+Co G Hi(f2_, A) andc/^Q = (^q +co^ G 

3 

Ha A). Moreover the choice of cq gives the compatibility condition (3.30) of problem 
(3.28)for A; = 1. 

3 

Again we take ipi = (pi + Ci & H2 (1]„, A), with ip^ uniquely determined by (3.28) 
with k = 1 under the condition (p^ dx = 0. Then iff = (pf + ciip G H2 {Q^, A), with 
(ff uniquely determined by (3.29) with (pf = on S and 

(3.34) ci = - J d„(pt (^s / J d„^ds. 

3 

We can continue this iterative process to construct the sequences {v^^T }a;>o C H2 (f2„, A) 
and{v?+}fc^o C Hi(n+,A). 

3.3.2. Proof of Theorem 2.1 in the Dirichlet case. The absolute convergence of the series 
Sfc>o P~^'^k ™ -'^^ (^±5 ^) obtained like in the Neumann case. The proof of the unique- 
ness of solutions to problem (2.1) in the Dirichlet case is also similar to the Neumann case. 

3.4. Complements. In this subsection we give some complementary results, sharper or 
more general than those of Theorem 2.1. 

3.4. 1 . Uniform estimates for ip — tpQ. As a consequence of the bounds (3.20), we have 

n=l 

and we deduce the following estimate between the solution ip E V oi the problem (2.1) 
and the solution lpq of the limit problem as p tends to infinity. 

Theorem 3.2. Let us assume that a+ ^ 0. There exist a constant po > independent 
of such that for all a_ G {z G 'C\\z\ ^ pok+|}> the unique solution ip E V of the 
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problem (2.1) with data (/, g) satisfying (2.2)-(2.3) converges in the piecewise H^/^ norm 
to the solution ipo of the limit problem as p tends to infinity, with the uniform estimate 

(3.35) ||^+ - <^^|||,f^^ + 11^" - ^0 ^ C'pJpr^(|a+|-^||/||o,n + ||^||o,e) 

with a constant C > 0, independent of a^, a„, /, and g. 

In this context, the above result gives sharper estimates than [8] where we find a charac- 
terization of limit solutions and strong convergence results for similar (and more general) 
problems. 

Likewise, an estimate of the remainder at any order is valid: 

K K 
n=0 n=0 

3.4.2. Uniform estimates when a_ and a+ play symmetric roles. The framework of The- 
orem 2.1 can be extended so that a_ and a+ play symmetric roles, and so that the contri- 
butions of the norms ||/^||o,n_ and ||/"'"||o,n+ are optimally taken into account in the right 
hand side of estimates. For this, we require the following assumptions 



(3.36a) / /dx= / ^ds = li V = Vn 

(3.36b) / /+dx= / /~dx= [ gds = if \/ = Vd and |a_| > |a+| 



(3.36c) no condition if = Vd and |a_| ^ |a+|. 

Proposition 3.3. Let us assume that a± ^ 0. There exist a constant po > such that for 
all couples (a_, a+) such that 

|a_| ^ Pok+I or |a+| ^ pok-l 

the problem (2.1) with data {f,g) satisfying (3.36) has a unique solution ip E V, which 
moreover is piecewise H^/^ and satisfies the uniform estimate 

,n nn\ II +11 I II -II ^ IIO'^- , ll/'*'llo,!^+ 

(3.37) + y ^CpA— — ■ — + — — — + ||^||o,e 

V I I I ^+ 1 

with a constant Cpf^ > 0, independent of a_, f, and g. 

Proof 1) Let us first prove estimate (3.37) in the Neumann case and when the modulus of 
p := a_(a+)~^ is large enough. After the change of data (/, g) g) as explained 

at the beginning of subsection 3.1, proving estimate (3.37) reduces to show 

(3.38) ytJin^ + llv^^P)!!!.^- ^ ^(-"o) (l-^^'l Wf'Wo,^^ + Wf^hn^ + MW^)^ 

instead of (3.23) (note that the new factor |p~^| in front of ||/~||o,q_ is equal to |a+|/|a_|). 

Estimate (3.38) is in fact a mere consequence of estimates (3.21) and (3.22), where we 
take advantage of the presence of the factor |p~M in front of the norm llv^rlL a of ^pV 
in (3.21). 
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2) Still in the Neumann case, but when p := a+(a_)~^ is large enough, the sequence of 
problems to be solved is now 

= in Q+ 

(3.39) <; d„^+ = g on S 

= on dfl 




(3.40) 

and for = 1, 2, 



A(/)q = f in Q_ 
= fo on S 



= Slf+ in fi+ 

(3.41) { d„^l = -5lg + d„^l^, on S 



(3.42) 




= on dfl. 

Ay9^ =0 in r2_ 
= on S 

The compatibility of the right hand sides of problems (3.39) and (3.41) in can be 
checked by arguments similar to those used in the case when a^{a^)"^ is large (§ 3.2.1). 
The estimate can be proved similarly. 

3) In the Dirichlet case, if |a_| >> |a+|, under assumption (3.36b), we see that in (3.33), 
we simply have 



'^ = ~J dnif^dsj j dn'ipds. 



Thus / does not influence (po, and we have estimates like in (3.21) with the factor \p ^\ 
in front of ||q37 11 a , . We deduce estimate (3.37) like in the Neumann case. 

4) Finally, in the Dirichlet case, if |a_| << |a+|, none of the elementary problems is of 
Neumann type. Hence no compatibility condition is required and we can prove estimate 
(3.37) as previously. □ 

The compatibility conditions (2.3) (and a fortiori (3.36a)-(3.36b)) are not necessary for 
the solvability of problem (2.1): For Neumann exterior boundary condition, the necessary 
and sufficient condition is 



/ /dx+(a+ — a_) gds 
Jn is 



0. 



It depends on coefficients a±. If we want to have the compatibility of the right hand side 
for any value of the coefficients a± we can either assume (2.3) or replace the coefficient in 
front of the integral by 1, defining the new problem 

(3.43) \/ipeV, / a+V^+- V^+dx+ / a_Vv?"- V^'dx = - / /^dx+ / g^ds. 
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Proposition 3.4. If we assume the compatibility conditions 



(3.44) - / /dx + / gds = Q if V = Vu and nothing if V = Vd , 



Jn Jt, 



then problem (3.43) is uniquely solvable if the modulus of p := a_{a^) ^ is large enough 
and its solution satisfies the uniform estimate 



Remark 3.2. We pay the weaker assumption on the data by a weaker estimate since we 
have the same estimate for solutions of problem (3.43) as for solutions of problem (2.1): In 
problem (2.1) the interface datum is (a+ — a-)g (the coefficient |a+ — a_| tends to infinity 
as IpI — oo) whereas the interface datum of problem (3.43) is g alone. 

Proof. The construction of the terms of the series expansion is similar as in the proof of 
Theorem 2.1. Now we have = in the Neumann case, and v?o = in the Dirichlet 
case. □ 

4. Proof of uniform estimates for Maxwell solutions at high 

conductivity 

We consider now the harmonic Maxwell system (1.2) at a fixed frequency cu satisfying 
Hypothesis 2.2. We are going to prove the following sequence of statements: 

Lemma 4.1. Under Hypothesis 2.2, there are constants ctq and Cq > such that if a ^ ctq 
any solution (E, H) G L^(f2)^ of problem (1.2) with boundary condition (2.8a) and data 
j e Ho(div, Q) satisfies the estimate 

(4.1) l|E||o,n ^ C'o||j||H(div,Q)- 

A similar statement holds for boundary conditions (2.8b) and data j G H(div, Q). 

This lemma is the key for the proof of Theorem 2.3 and is going to be proved in the next 
subsection, using in particular our uniform estimates in the scalar case (this is the main 
difference with the proof of Theorem 2.1 in [7]). As a consequence of this lemma, we will 
obtain estimates (2.1 1): 

Corollary 4.2. Let a > 0. Let{E, H) G L'^{^)^ be solution of problem (1.2) with boundary 
condition (2.8a) and data] G Ho(div, Q). IfE satisfies estimate (4.1), then setting 



there holds 

(4.2) ||E||o,f7 + II curlE||o,n + || div{iujeo - a)E\\o^a + y/a \\E\\o^q_ ^ Ci||j||H(div,Q)- 
A similar estimate holds for boundary conditions (2.8b) and data j G H(div, Q). 
Finally, estimate (4.1) implies existence and uniqueness of solutions. 

Corollary 4.3. Let a > 0. We assume that estimate (4.1) holds for any solution (E, H) G 
L^(r2)^ of problem (1. 2)-(2.Sa) with] G Ho(div, f2). Then for any ] G Ho(div, fi), there 
exists a unique solution (E, H) G L^(r2)^ of problem (1.2)-(2.8a). A similar result holds 
for boundary conditions (2.8b) and data] G H(div, Vt). 




+ \\v-\\ln^^CJa4-\\\f\\ 



o,f7 + ||fi'||o,s) • 
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The previous three statements clearly imply Theorem 2.3. To prepare for their proofs, we 
recall variational formulations in electric field for the Maxwell problem (1.2) with bound- 
ary condition (2.8b) or (2.8a), cf. [13] for instance. Let 

(4.3a) H(curl, fi) = {u G \ curl u G L2(fi)} 

(4.3b) Ho(curl, fi) = {u G H(curl, Q) \ u x n = on dQ}. 

If (E, H) G L2(fi)2 is solution of (1.2)-(2.8a), then E G H(curl,l]) satisfies for all E' G 
H(curl,fi): 

(4.4) / ( curl E- curl E^ - K^E ■ E^) dx-zi^a / E ■ ^ dx = zz/ / j ■ E^dx 

where we have set k = cuy/eofxo and = w/iQ. If boundary conditions (2.8b) are consid- 
ered, then E G Ho(curl, fi) and (4.4) holds for any E' G Ho(curl, fi). 



4.1. Proof of Lemma 4.1 : Uniform L^ estimate of the electric field. Reductio ad ab- 
surdum: We assume that there is a sequence (E^, H„i) G L^(r2)^, m G N, of solutions of 
the Maxwell system (1.2)-(2.8a) associated with a conductivity am and a right hand side 
G Ho(div,fi): 



(4.5a) curl E^ — iw/ioHm = in ^2 , 

(4.5b) curl + {iueo - a^) E^ = j„ in n, 

(4.5c) H„, X n = on dQ , 

satisfying the following conditions 

(4.6a) (Tm — 7- oo as m — 7- oo, 

(4.6b) ||E„J|o,n = l VmGN, 

(4.6c) ||jmllH(div,n) -^0 as m ^ oo. 



Note that the external boundary condition E^ ■ n = on dQ is but a consequence of the 
equation (4.5b), the boundary condition (4.5c) and the condition j ■ n = on dil contained 
in the assumption that belongs to Ho(div, f2). 

We particularize the electric variational formulation (4.4) for the sequence {Em}: For 
all E' G H(curl,fi): 

(4.7) / ( curl Em ■ curl E' — k^E^ • E') dx — iuam / E^ ■ E' dx = zi^ / j^, ■ E' dx . 

Jn Jn^ Jn 

Choosing E' = Em in (4.7) and taking the real part, we obtain with the help of condition 
(4.6b) the following uniform bound on the curls 



(4.8) II curlEmllcn < K,+ V^^lllmllo,!^ • 
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4.1.1. Decomposition of the electric field and bound in H^. We recall that we have as- 
sumed that the domain Vt is simply connected and has a smooth connected boundary. Re- 
lying to Theorem 2.9 and Theorem 3.12 in [2], we obtain that for all n G N there exists a 
unique w,„ e H^(r2)'^ such that 

(4.9) curlWm = curl E^, div = in fi, and ■ n = on dVt . 
Moreover, we have the estimate 

(4.10) ||w„,||i,j7 ^ C|| curlE^llcn, 

where C is independent of m. As a consequence of the equality curl Wm = curl Em and 
the simple connectedness of f2, we obtain that there exists y?^ G H^(r2) such that 

(4.11) Ern = Wm + Vipm- 

We write equation (4.5b) as 

curl Hm + {iOJEo - £m)(Wm + Vv?m) = im- 

Let ip E H^(n) be a test function. Multiplying the above equality by Vijj and integrating 
over f2, we obtain, using that div = 0: 

(4.12) / {iuEo - a^) Vipm -V^p dx = - / divj^ tp dx - / w„ ■ n |^ ip ds . 

Note that the boundary values curl Hm ■ n = and ■ n = on dVt have been used here. 

Thus ifm is solution of the Neumann problem defined by the variational equation (4.12). 
Since 

divj^ e L^(ri) and / divj„dx= / ■ n ds = 0, 

Jo. Jan 

and 

Wm ■ n G L^(S) and / w„i ■ n |^ ds = / divw^dx = 0, 

the Neumann problem defined by (4.12) satisfies the assumptions of Theorem 2.1 with 
a_ = itueo — cr,„ and a+ = iusQ. Therefore we have the following uniform estimate for 
am large enough (i.e. for m large enough, cf. (4.6a)) 

lkmll|,n+ + Ibmllf,f7_ ^ C'odi divj„||o,f7 + ||w„ ■ n||o,s). 
Since ||wm ■ n||o,s is bounded by ||wm||i^n, the above inequality implies 
(4-13) + Wv^WlM- ^ "^odl divj„||o,n + ||w„,||i,f7). 

Finally (4.6c), (4.8), (4.10) and (4.13) implies that 
(4-14) \Wtr\\ln+ + ll^mllf,n_ + \\^Ai,n ^ B 

for a constant i? > independent of m. With (4.1 1), (4.14) gives that the sequence {Em} 
is bounded in H2 on and 

l|E™lli,n+ + l|Emlli,n_ ^B. 
Combining the above bound with (4.8), we obtain the uniform bound 
(4.15) l|E+||i n+ + l|E- 111 n_ + II curlE„||o,n ^ C. 



(4.16) 
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4. 1 .2. Limit of the sequence and conclusion. The domains VL± being bounded, the embed- 
ding of H2 in L^(f2-|-) is compact. Hence as a consequence of (4.15), we can extract 
a subsequence of (still denoted by {£„}) which is converging in L^(fi). By the 

Banach-Alaoglu theorem, we can assume that the sequence curl is weakly converging 
in 'l?{Vt): We deduce that there is E G L^(l]) such that 

curl Em ^ curl E in L^(r2) 
Em ^ E in 

A consequence of the strong convergence in L^(f2) and (4.6b) is that ||E||o,f7 = 1. Using 
Hypothesis (2.2), we are going to prove that E = 0, which will contradict ||E||o,c = 1, and 
finally prove estimate (4.1). 

Taking imaginary parts in (4.7) when E^ is the test-function, then letting m — )■ +00 and 
using (4.6c) we get ||E||o,n_ = 0. Hence, 

(4.17) E = in Vt_ . 

Let us introduce the space 

Ho(curl, S) := {u G H(curl, fi+) | u x n = on E}. 

In particular, (4.17) implies that E"*" := E belongs to Ho(curl, S). 

Let $ G Ho(curl, E). Then the extension $0 of $ by on defines an element of 
H(curl, Vl). We can use $0 as test function in (4.7) and we obtain 

/ (curlEm ■ curl$ - K^Em ■ $) dx = iz/ / j^-^dx. 

According to (4.16) and (4.6c), taking limits as m — +00, we deduce from the previous 
equalities 



(4.18) 



/ ( curl E+ ■ curl $ - k'^E+ • $) dx = 
Jn+ 



i.e., E+ G Ho(curl,fi+,E) satisfies (4.18) for all $ G Ho(curl, E). Integrating by 
parts we find (with $t the tangential part of $ on d^l) 

(curl E+,curl$)o,n+ = (curl curl E+, <l>)o,f7+ - (curl E+ x n, <l>T)|aQ. 

Thus we have 

{curl curl E"^ — E'^ =0 in 
E+ X n = on E 

curl E"*" X n =0 on dfl. 

Setting := (iw/io)"^ curl E"^, we obtain that curl H"^ = —iujeoE~^ and we deduce the 
remaining boundary conditions 

H+ ■ n = on E and E+ ■ n = on 9fi 

from the previous relations. Hence (E'^, H'^) G L^(r2+)^ is solution of problem (2.9). By 
Hypothesis 2.2, we deduce 

E+ = in Q,. 
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Hence, with (4.17), we have E = in fi, which contradicts || E||o,c = 1 and ends the proof 
of Lemma 4.1. 

4.2. Proof of Corollary 4.2. Let (E, H) G L'^{^1)'^ be a solution of the Maxwell problem 
(1.2) with boundary condition (2.8a) and data j G Ho(div, ^l). We assume that 

(4.20) l|E||o,n ^ C'o||j||H(div,f7)- 

Then E G H(curl, f2) is solution of the variational problem (4.4). Taking as test function 
E itself, we obtain the identity 

(4.21) / ( curl E- curl E - K^E ■ E) dx-iz/cr / E ■ E dx = iz/ / j ■ E dx . 

Taking the real part of (4.21), we obtain 

II curlE||2^f^ = K^||E||^^f^ - u lm(j, E)o,n 
hence, using inequality (4.20) and Cauchy-Schwarz inequality. 



(4.22) II curl E||o,f, ^ {kCo + ||j||H(div,f7)- 

Then, taking the imaginary part of (4.21), 



c^l|E||oo_ = -Re(j,E)o,: 



hence, 

(4.23) v^||E||o,n_ ^ llj||H(div,f^)• 
Taking the divergence of equation curl H + {iueo — ct)E = j, we immediately obtain 

(4.24) II dw{iujeo - o:)E||o,f7 = || divj||o,n. 
Formulas (4.22) to (4.24) yield Corollary 4.2. 

4.3. Proof of Corollary 4.3. Let a and tu (i.e., k) be fixed. Let us introduce the piecewise 
constant function aonQ 

(4.25) a = l + —a . 

With this notation, the sesquilinear form in the left hand side of (4.4) becomes 

(4.26) I ( curl E ■ curl £! - K^aE-E!) dx. 

The proof of Corollary 4.3 relies on a classical regularization procedure: We consider the 
functional space 

XT(a) = {E G H(curl,fi)| div(aE) G 1^(1]), E ■ n = on 9^}. 

Let s > be a real number, which will be chosen later. Let us introduce the sesquilinear 
forms As and B : y^j{q) x Xt(«) — ^ C 

(4.27a) A,(E, E') = / (curl E ■ curl E^+ s divaE div^) dx 

Jo. 

(4.27b) 5(E,E') = / aE- E^dx. 
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With a right hand side j G Ho(div, fi), we associate a new right hand side depending on 
the parameter s defined as an element of XT(a)' by 

(4.28) j,(E') = / (j ■ r - 4 divj div^) dx VE' G Xjia). 

The regularized variational formulation is: Find E G XT(a) such that 

(4.29) VE'gXt(«), ME,E')-K^B{E,E') = tul{E'). 
As a consequence of [4, Th. 7.2], we obtain that if 

(4.30) — is not an eigenvalue of the Neumann problem for the operator div a V, 

s 

then any solution (E, H) G L^(r2)^ of problem (1.2)-(2.8a) with j G Ho(div, fi) provides a 
solution of problem (4.29), and conversely, any solution E of (4.29) provides a solution of 
(1.2)-(2.8a) by setting H = (iw/i)-^ curl E. 
Thus, we choose s so that (4.30) holds. 

Since the form Ag is coercive on Xj{a) and the embedding of L^(f2) in XT(a) is com- 
pact, we obtain that the Fredholm alternative is valid: If the kernel of the adjoint problem 
to (4.29) 

(4.31) Find E' G Xt(«), VEGXjla), v4,(E, E') - k^5(E, E') = 0, 

is reduced to {0}, then problem (4.29) is solvable. 

We see that the assumption of Corollary 4.3 implies that (4.3 1) has only the zero solution, 
and that the same holds for the direct problem 

(4.32) Find E G Xt(«), V E' G XT(a), A,(E, E') - k^B{E, E') = 0, 
of course. 

All this implies the unique solvability of problem (1.2)-(2.8a) with j G Ho(div, Vl). 

5. Application: Convergence of asymptotic expansion at high 

conductivity 

In the Maxwell case, see equations (1.2), let us introduce the parameter 

Thus, when a — )• oo, 5 tends to 0. Note that the function a defined in (4.25) can be written 
(5.2) a = ln+ + (l + ^yn_. 

Several works are devoted to the interesting question of an asymptotic expansion as 5 — )• 
of solutions of the Maxwell system (1.2) with complementing boundary conditions on d^l 
when the interface S is smooth: See [15, 10, 11] for plane interface and eddy current 
approximation, [7] for impedance boundary conditions and [14] for perfectly insulating or 
perfectly conducting boundary conditions. 
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5.1. Assumptions. We assume that S is a smooth surface, and we follow the approach 
of [14]. In order to fix ideas, we take perfectly insulating boundary condition (2.8a) and 
assume Hypothesis 2.2 for this condition. By Theorem 2.3 there exists such that the 
conclusions of the theorem hold. From now on we assume that 



(5.3) 



cr ^ (Jo, i.e. 5 ^ §0 with Sq 



UJEq 



Let j G Ho(div, fi) such that j = in Then for all 5 ^ 5o, there exists a unique 
solution to problem (1.2)-(2.8a), which we denote by (E(5), H(5)). Then it is possible to 
construct series expansions in powers of 6 for the electric field Ej^ in the dielectric part 
fi-i. and E^-) in the conducting part 



(5.4a) 
(5.4b) 



J>0 



with E7(x;5) = x(y3)Wj(2/^, 



In (5.4b), y = (y^, y^) are "normal coordinates" to the surface S in a tubular neighborhood 
U- of E in the conductor part In particular, represents the distance to S. The 
function y i-)- xiva) is a smooth cut-off with support in U_ and equal to 1 in a smaller 
tubular neighborhood of S. The functions \Nj are profiles defined on S x IR+. Moreover, 
for any j eN 



(5.4c) 



E+ e H(curl, n+) and W^- e H(curl, E x R+) 



There hold a similar series expansions in powers of 5 for the magnetic field ^(S)- 

The validation of the asymptotic expansion (5.4) consist in proving estimates for remain- 
ders Rm- 5 defined as 



(5.5) 



■(5) 



3=0 



in . 



This is done by an evaluation of the right hand side when the Maxwell operator is applied 
to Rm;5. By construction [14, Proposition 7.4], we obtain 



(5.6) 



m; i5 



curl curl - ti^a+R^.^ 
curl curl R".^ — K^a_R; 

[ curl R„; 5 X n] ^ 
curlR+.5 X n 



= 


in 




Jm; 5 


in 




= 


on 


E 


8m; <5 


on 


E 


= 


on 


dQ 



Here, according to (5.2), a+ = 1 and a_ = 1 - 
E X n across E. The right hand sides (residues) j^. g 



and [E x n]2 denotes the jump of 



and g^. g are, roughly, of the order 5^ 
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5.2. Convergence result. The main result of this section is the following. 

Theorem 5.1. Under Hypothesis 2.2 in the framework above (section 5.1), we assume that 
we have for all m the following estimates for the residues j~. g and g^. g in (5.6) 



(5.7) \\'im:shn- + llgn^;^lli,s + l|curlsg„.5||3 s ^ 



(5.10) 



where Cm > is independent of 6, and rrio G N independent of m and S. Then for all 
m E N and 6 G (0, 6o], the remainders Rm:S (5.5) satisfy the optimal estimates 

(5.8) \\Ri,shn+ + ||curlR+.5||o,n+ +5"^||R~;5llo,n_ +5^||curlR-.5||o,n_ ^ C'J'^+K 

Proof Here we denote by C,„ various constants which may depend on m but not on 5. 
Step 1. We cannot use Theorem 2.3 directly because curl curl Rm-s — «^^ttR»n;(5 does not 
define an element of H(div, fi). We are going to introduce two correctors C^; s and s 
satisfying suitable estimates and so that 

{[{R.m;5 - Cm;s) X =0 ou S 

[curl(R„;5 - C^;5) X n]^ =0 on S 

curl(R^.5 - C„;5) X n =0 on dQ 

and 

[«(Rm;5 - C,„;5 - D„,.5) • n]^ =0 on S 
[(R^;5 - C^;5 - D^.^) X n]^ =0 on S 
[curl(R^.5 - C„;5 - D„.5) X n]^ =0 on S 
^ curl(R„.5 - C„;5 - 0^.5) X n =0 on 9fi . 

Step 1a. Construction of C„i;5'- We take Cm;S = in and use a trace lifting to define 
Cm- 5 in It suffices that 

' C;5 X " =0 on E 

(5.11) { curlC+.^xn = g^.^ on E 

curl C^. g X n = on dQ . 

Denoting by and C3 the tangential and normal components of C,"^.^ associated with a 
system of normal coordinates y = (y^, y^), and by g/^ the components of g^.^ the above 
system becomes (cf. [14, Proposition 3.26]) 

f = on E 

(5.12) I dsC^-d^Cs = gp on E 

y d^iCfi-dfiCz =0 on dn. 

It can be solved in H^(f2+) choosing C3 = and a standard lifting of the first two traces 
on E and dVt with the estimate 

(5.13) l|C+ ,||2,n+ ^CIIg^^^lliE- 
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Step 1b. Construction of D^.^: Let us denote Rm;5 — ^m;5 by S for short. Again, we 
take D".^ = and use a trace lifting to define D^.^. It suffices that 

f 0+5- n = [aS-n]^ on S 

(5.14) <^ D+^xn =0 on S 

[curlD+.^xn =0 on S U 91] 

In normal coordinates and associated components, these conditions become, compare with 
(5.12) 

{Ds = [aS ■ n] ^ on S 

Dp =0 on S 

d-iDp-dpD^ =0 on 

which can be solved in H^(f2+) (first determine D^, then Dp) with the estimate 

(5.16) ||D+,||2,f,^ ^C||[«S-n]^||3,s. 

Since [ curl S x n] ^ = 0, we find that, by construction 

—K^ \aS -11]^= [(curl curl S — h?qS) ■ n] ^ 

= - (curl curl C+. 5 - k^C^.s) Is ' " " V^ls ' " 

= curls - j;;.^!^. ■ n . 

Hence 

(5-lV) ||[aS-n]^|||_5. ^ ||curlsg^.5|||_s+ W'hn.sh^- ■ 

We deduce from assumption (5.7), and (5.13), (5.16), (5.17) 

l|C+,i|2,c+ + ||D+,|h,n, 

Since by construction C~.^ = D^.^ = and \Crn;S ^ "]s ^ [D^.^ x n] ^ = 0, the above 
estimate implies 

(5.18) ||C„,5||o,n + II curlC„;5||o,f7 + ||D„,;5||o,f7 + || curl D^.^llo,^ ^ C^^""-™" . 
We set 

(5.19a) Rm;5 := Rm;5 — Cm;5 — Dm;5 

and 

(5 . 1 9b) '3m;S-= curl curl R„. s - K^a^m; s ■ 

Hence by construction, ^ G H(div, f2) with the estimates 

(5-20) |bm;5llH(div,Q) ^ CmS''' 



irm—rno 



Step 2. We can apply Theorem 2.3 to the couple (E, H) = (Rm;5, ^ curl Rm^s) 

and, thanks to (5.19b), obtain 

||Rm;5||o,n + II curl Rmj^llo.Q ^ C*|| Jm; 5 II H(div,f7) • 
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Combined with (5.20), this gives 

(5.21) ||R„^;5||o,^ + II curlR^;,||o,Q ^ C^6"\ 
Together with (5.18) and (5.19a), this estimate gives finally 

(5.22) \\R^;5\Wn + II curl R^.,s\Wu ^ C^5"\ 

Step 3. In order to deduce the optimal estimate (5.8) for R^; 5, we use (5.22) for m + 1 + 
mo, which yields 

(5.23) II Rm+l+mo;(5||o,Q + || CUll Rm+l+mo; 5 1| 0,^ ^ Cm^^^^- 

But we have the formula 

m+l+mo 

(5.24) Rm;5 = ^ ^■'Ej + Rm+i+mo;5- 

j=m+l 

Moreover by definition the E^^ do not depend on 5 and the are profiles: Using (5.4b) 
we find that for any j G N 

(5.25) r2||E-||o,f7_ +5^||curlE7||o,Q_ ^ C||W,|| H(curl,SxR+)- 
Combining with (5.4c), we obtain for any j G N 

(5.26) IIE+II 

H(curl,n+) 

We finally deduce the wanted estimate (5.8) from (5.23) to (5.26). □ 

Remark 5.1. As a consequence of the works [7] and [14], we find the existence of asymp- 
totics of the form (5.4) when the interface S is smooth, if the right hand side j is smooth 
and has its support in the dielectric part Moreover, estimate (5.7) is true for mo = 1, 
cf. [14, Ch.7]. 

Remark 5.2. If the interface has conical points, or is polyhedral, many difficulties are en- 
countered for an asymptotic analysis. We refer to [12] for an investigation of a scalar 
transmission problem with high contrast in polygonal domain. 
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